Quantum oscillations in the mixed state of d-wave superconductors 



O 

o 

(N 

> 

O 



O 

o 

I 

Oh! 



I 

c 
o 
o 



> 

o 



o 



X 



Ashot Mclikyan 1 and Oskar Vafck 2 

1 Materials Science Division, Argonne National Laboratory, Argonne, IL, 60439 

2 National High Magnetic Field Laboratory and Department of Physics, 

Florida State University, Tallahassee, Florida 32306, USA 

(Dated: February 2, 2008) 

We show that the low-energy density of quasiparticle states in the mixed state of ultra-clean 
d— wave superconductors is characterized by pronounced quantum oscillations in the regime where 
the cyclotron frequency ftw c <C Ao, the d— wave pairing gap. Such oscillations as a function of 
magnetic field B are argued to be due to the internodal scattering of the d-wave quasiparticles near 
wavevectors (doko, ±&d) by the vortex lattice as well as their Zeeman coupling. The periodicity of 
the oscillations is set by the condition ko\/hc/(eB) = Icd' \/hc/(eB') (mod 2ir). We find that there 
is additional structure within each period which grows in complexity as the Dirac node anisotropy 
increases. 



The behavior of <i-wave superconductors in magnetic 
field has recently come into sharp focus due to the exper- 
imental observation of the oscillations in the longitudinal 
and Hall electrical transport at extremely high magnetic 
fields for very clean underdoped YBa2Cu3C>6.5 [l| and 
YBa 2 Cu40g [2|, y)]. At temperatures about 2 K and for 
magnetic fields in excess of about 30 T, the Hall resis- 
tance is finite, negative (i.e. electron like), and oscillates 
as a function of magnetic field, with approximately four 
pronounced peaks between ~ 45 T and 60 T. Below about 
30 T, the system is in the mixed superconducting state, 
and the electrical Hall conductivity vanishes. The cur- 
rent interpretation of these oscillations involves charge 
or spin order induced reconstruction of the Fermi surface 
with electron pockets ( [l|, y, |j , [5j, [6J] ) , which result from 
the appropriate tuning of the strength of the charge/spin 
potentials in the particle-hole channel [4j, [5|, [6( . 

As the estimated mean-field H c2 for the YBCO sam- 
ples in question significantly exceeds 60 T[7j, it is at 
present unclear why one should expect the superconduct- 
ing order parameter amplitude to collapse at the fields 
applied in Rcfs. [l|, |2|, |3|] . Instead, the local superconduct- 
ing correlations should persist, and rather the system is 
expected to be in the vortex liquid state [3, |9(. Addi- 
tionally, the sign reversal of the Hall effect, from positive 
to negative, with decreasing temperature near T c in op- 
timally doped YBCO has been observed and interpreted 
earlier as resulting from the flux flow|l0l. I 111 Il2| . The 
structure factor of such liquid state is expected to closely 
resemble that of a vortex solid. In this context, the nat- 
ural question is the existence, the periodicity, and the 
physical nature of such magnetic field induced oscilla- 
tions in the vortex state. 

In this work we therefore examine the properties of 
the d— wave quasiparticles (qps) in the vortex solid state 
at intermediate magnetic fields, i.e. we examine the 
particle-particle channel. For the square vortex lattice 
observed in the small angle neutron scattering [l3j, we 
find that the combined effect of the orbital and spin cou- 
pling of the qps to the magnetic field does induce os- 



cillations in the low energy density of states, N(E), and 
that the nature of these oscillations depends on the Dirac 
cone anisotropy ao- In the physically relevant regime 
where the d-wave gap Ao 3> h\o c (the cyclotron frequency 
oj c = eB/mc), the density of states N(E) is an oscillatory 
function of kp£, where ko is the k— space half-distance 
between the nearest nodes (see Fig. 1) and £ is the mag- 
netic length £ = ^Jhc/ Be. For ao = 1, the period of 
the oscillations corresponds to 8(kr)£) — 2n, but with 
increasing anisotropy the number of oscillations within 
the period increases. In the absence of Zeeman cou- 
pling N(E) is exactly periodic in the scaling regime and 
the spectrum oscillates between fully gapped and nodal. 
Upon inclusion of the Zeeman coupling, which cannot be 
neglected for fields of Ref. [l|, |2j, y] , the oscillatory be- 
havior persists, but the scaling is only approximate. For 
au = 1 and for the fields in the excess of 30T, the com- 
bined effect of the orbital and Zeeman coupling leads to a 
sequence of transitions and the system oscillates between 
a thermal metal with finite N(0) and a thermal insulator 
with vanishing N(0). A similar pattern of oscillating be- 
havior holds for large an , albeit with effectively increased 
frequency and significantly decreased gap. The envelope 
of the N(0) oscillations follows an approximately £~ x be- 
havior. These findings should be contrasted with the os- 
cillations found in the regime of extreme high fields Ref. 
14j where the number of Landau levels below the Fermi 
energy range from 1 to 7. 

The starting point is the lattice Bogoliubov-de Gennes 



(BdG) eigenequation[15l. Il 



WoV'r = E*I>T 



(1) 



where the Hamiltonian TLq acts on the two component 
Nambu spinor ip r = [u r ,z» r ] T and has the following ex- 
plicit form 



TLa 



£ T — n + ^hbB 

a: 



-£* + fi + %/j.bB 



(2) 



Both £ r and A r are defined through their action on a 
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The physical supcrfluid velocity enters via the factor 



FIG. 1: Left: Magnetic unit cell £ x £ of a tight-binging lattice 
(solid lines) containing two vortices joined by a branch-cut 
shown for £ = 6a. Right: Fermi surface (black solid lines) 
of cuprate superconductors with four nodal points at k^! = 
(±fc_D, ±fco) in a Brillouin zone of size 2-n/a. In the presence 
of the vortex lattice, the Brillouin zone is reduced to 2ii /£ 
(spacing of the blue grid). The properties of the low energy 
states depend on the commensuration between the inter-nodal 
distance 2fcjv and the reduced reciprocal lattice vector 2n/£. 



lattice function f r as 



t 



-iA r 



' fr+S, 



'Vsfr+S- 



(3) 
(4) 



<5=±x,±y 

A r / r = A Y, e 

i5=±x,±y 

In the symmetric gauge, the magnetic flux $ through 
an elementary plaquctte enters the Peierls factor via 
A rr +£ = — 7ry$/0 o > A rr+ y = 7rx$/0o; the electronic flux 
quantum is 0o = hc/e. The d-wave symmetry is encoded 
in rjs = +(— ) if 5 || x(y). Importantly, 8 rr i winds by 2tt 
around each of the magnetic field induced vortices. The 
initial Ansatz[17l] for the pair phases is 

e 10 ™' = (e* + e # - 



'VI 



,i<i>* 



+ e 



i(i) r 



(5) 



where V x V</>(r) = 2ttz £\ S( T - r *) and v ' W(r) = 
where r^ denotes the vortex positions. Although the 
phase of the order parameter does depend on the pre- 
cise form of the self-consistency condition, the deviations 
from the adopted form are weak, and moreover, both the 
symmetry of the phase and its singular part are fixed 
unambiguously by the vortex lattice. 



Connecting pairs of vortices by branch cuts 17 1 , we 
can define the singular gauge transformation 17J, [l8| 
U = e^' T3 ^ r where the Pauli sigma matrices act on the 
Nambu spinors. The transformed Hamiltonian 7i(k) = 



-ik-r 



U- 1 H Ue ikr becomes 
W(k) = a 3 (&(k) - A') + oiA r (k) + Z^bB, 
where the transformed lattice operators satisfy 



(0) 



£ r (k)^ r = 


= -t J2 Z2,rr+sxe i ° sV "+ s e ik - s i> r+s (7) 




<5=±x,±y 


A r (k)^r = 


= A Y] z 2 , rr +s x r)se tk ' s i) r+s . (8) 



iV„, _ „-iA„, 



e " 



e~ iA "' (1 + e 1 ^'-^)/]! + e^'-^l (9) 



and represents the lattice analog of the semiclassical 
(Doppler) effect [19l|. The Z2 field Z2,rr' = 1 on each 
bond except the ones crossing the branch cut where 

Z2,rr> = -1. 

In the magnetic fields of interest the vortices of YBCO 
form a square lattice with primitive vectors oriented 
along the d-wave nodes [13J, and therefore the trans- 
formed Hamiltonian ^ is invariant under discrete trans- 
lations by the primitive vectors Ri = £x and R2 = ly 
defining the magnetic unit cell, reflecting the periodicity 
of Vrr' and the periodic choice of the branch cuts. Con- 
sequently, it can be diagonalizcd in the Bloch basis. By 
the Bloch condition 7i(k) acts on the periodic functions, 
and the crystal wavevector k varies continuously within 
the I s * Brillouin zone defined by the primitive reciprocal 
lattice vectors Ki = 2ttj, K 2 = 27r|. Note that, since 
the unitary transformation hi is time independent, the 
Hamiltonians |(5J) and (J5|) have the same thermodynamic 
and tunneling density of states. 

The above claims follow from our detailed numerical 
study of the tight-binding Bogoliubov-de Gennes Hamil- 
tonian ^. Since the Zceman term is just a simple overall 
shift of the qp energies, let us first neglect it. The results 
of the numerical diagonalization of ^ can be summa- 
rized by the following scaling form of the qp eigenener- 
gies: 



E n (k) 



hvt 



T n (k£,a Dl k D £) 



(10) 



<5=±x,±y 



Here vf and va are the Fermi and gap velocities of the 
d— wave qps near the nodes (±A;£>,±fc.D) (i n our m °dcl 
dD) vf = 2\/2athr 1 sin(fci)a) and Va = VpAo/t); T n is 
a dimcnsionless scaling function, which differs from the 
Simon-Lee scaling function in the following important 
aspect: it depends on krj£. Specifically, in the scaling 
limit of k]j£ — > 00, the function J- n does not approach a 
uniform value. This is amply illustrated in the Fig.([2|), 
where we used a color density plot to represent the den- 
sity of states per area rescaled as vf^N{E) for different 
values of ku(./{2-K) and E£/(Hvf)- 

As best seen for an = 1, near E = the spectrum 
is gapped although for krji w ir(2n + 1) the gap is very 
small, vanishing only at discrete set of points as expected 
for the unitary class |L6||. Moreover, the gap scales as £~ l 
and oscillates as a function of fcjj£ Physically, the oscil- 
lations are due to the strong internodal scattering which 
arises from the commensurability of the d-wave nodes 
and the vortex lattice. As seen from the three differ- 
ent panels, the pattern of the oscillations depends on the 
anisotropy an- With increasing arj, the DOS acquires 
ever richer structure with multiple minima and maxima 
- for presentation purposes we show only three periods 
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FIG. 2: (a-c): Oscillations of the qp density of states N(E) rescaled as lvFN(£E/(hvF)) are shown as a color map for old = 1 
(a), od — 4 (b), and an — 7 (c). The green line represents the shift of the qp zero energy due to the Zeeman term for t = 380 
meV, a — 3.8 A, g = 2, ko — 0.387r/a, and £ = 26a. (e-f): Zero energy density of states of ((2)1 as a function of magnetic field 
for a a — 1 (e), 4 (f), 7 (g) deduced from the figures of the top row and the generalized scaling form of eigenvalues (|10|l is shown 
in red dashed line. The solid black line corresponds to an effective broadening due to the inter-layer coupling 2t z cos(fc z c) with 
t z = 0.05t. 



(separated by white vertical lines) for ao = 4 and two 
periods for ao = 7. Additionally, as <xd increases, the 
scaling limit Ao 3> fno c is reached for smaller magnetic 
fields (larger €) and the approximate low energy scaling 
(fTTj)) holds for magnetic length ija ^$> 2-k^/od- 

The magnetic fields of experimental interest (45—70 T) 
correspond to £ € (20a, 25a), which is the regime where 
the Zeeman term cannot be neglected. Nevertheless, it is 
easy to take it into account, since it corresponds to a sim- 
ple shift in zero of the qp energy, represented by the green 
lines in each of the color DOS maps in Fig. [2] The DOS 
along the Zeeman "slice" is shown in Figs. [UJe-f). For 
anisotropy ar> = 1 , the spectrum exhibits pronounced os- 
cillations as a function of magnetic field (or kr>), changing 
from gapped (thermal insulator) with activated temper- 
ature dependence of the specific heat to gapless (thermal 
metal) with T— linear specific heat. This S£ = 2-KJkr> 
periodic sequence of transitions between thermal metal 
and thermal insulator is significantly richer for ar> > I . 
Within each such "primitive" period, a complex struc- 
ture develops, and the number of minima and maxima 
depends on ao- For ao = 7 and kn = 0.387r/a, which 
are approximately [2TJ the physical values expected for 
the YBCO samples of Ref. [1(, there are four "oscilla- 
tions" in the field range ~ 45 — 60T. Remarkably, this is 
precisely the experimentally observed " periodicity" [l| . 



Note that the effect described here occurs deep in 
vortex state with a fully developed superconducting 
amplitude and does not rely on scmiclassical orbits 
around electron pockets. Rather, as we now argue, 
it is a consequence of the commensuration involving 
the vortex lattice and the internodal separation (see 
Fig. [l}. In the scaling regime of interest, Ao ^ hw c , 
the appropriate starting point should be the linearized 
approximation 18l |21| of BdG Hamiltonian ©. In this 
approach the low-energy qps with E <C Ao are inde- 
pendent and behave as massless Dirac fermions inter- 
acting with vortices via the Doppler shift originating 
from the superflow v(r) = §V</> — -A. Additionally, 
the wavefunctions contain branch-cuts connecting vor- 
tices pairwiseJ3JT3JllJp, IHIH I3- The effective 
Hamiltonian [llliiniinijis then 



Ti-Un/VF 



Py 



V2 



03 + OD 



Py 



v/2 



Px V; 

o"i +m— 



+ v n 



V2 



(11) 



In this approximation, Tiu n exhibits the Simon-Lee [2 1| 
scaling Hu n (r,v F ,v A ) = ^Hu n (r/e,a D , 1) and the qp 

spectrum has a scaling form -E„(k) = ^p-.Tn (k^, Od). 
Now, we consider the corrections due to the terms left 
out during linearization. Provided that r is not in the 
immediate vicinity of a vortex core (r > £), the leading 
non-linear corrections, such as mv 2 (r)/2, are typically 



omitted on the basis that they arc smaller than the terms 
retained in (fTTj) by a factor of (fco-^) -1 which is small for 
typical fields of interest. However, such estimate is incor- 
rect. The effect of the non- linear terms is amplified [24 1 
by an anomalously large (low energy) qpwavefunctions, 
growing as r -1 / 2 near vortex locations [16J, [25|. This di- 
vergence is eventually cut off only at r « £. A typical 
"perturbation", such as mv 2 (r), results in the following 
matrix element between two eigenstates of TLu n at nodes 



,(j) 



(./') 



j and /: I = (*nkl mv ^ r )l*n'W- Due to Bloch symme 
try of the wavefunctions ^, the Bloch momenta k 



U) 



,0") 



and k^ ; +k' must differ by a reciprocal lattice vector G. 
Importantly, |/| oc ^r 1 ^MCi + C 2 cos(R ■ G)], which is 
of the same order as the terms retained in (|11[) . Here R 
is the primitive vortex lattice vector and the vortex core 
size £ serves as the small-distance cut-off of the other- 
wise divergent integrals [24j. Coefficients C\ t 2, in general 

depend on k and k', but not on k^ or k.p '. 

Thus, due to the r- -1 / 2 -behavior of the low-energy 
wavefunctions near vortex cores, the matrix elements I, 
and the resulting corrections to the energies, scale with 
magnetic length precisely as the energies of (|f 1|) , namely 
as ~ l~ l . The relative magnitude of the corrections 
due to the non-linearities compared to (jlip is determined 
by the magnitude of the parameter (fco^) -1 rather than 
(fco£) -1 . In cuprate superconductors the former is typi- 
cally of 0(1) since £ is of the order of a few lattice spac- 
ings. Thus, if the parameters of the BdG Hamiltonian (J6j) 
such /jorl are varied, the resulting qp spectrum evolves 
in a manner not captured by pip . 

Based on the above argument, for the square vortex 
lattice under consideration we expect the spectra for 
nodal momentum kjj (kr> ) and magnetic length i (£') 
to be similar when 



k D l = k D '£' (mod 2tt) 



(12) 



Therefore, the Simon-Lee scaling should be generalized 
to Eq. (flQ|) . In the regime where the scaling (fj"0|) holds 
the oscillatory part of the dispersion is fully determined 
by the product krji, and therefore the dependence of the 
spectrum on the magnetic field can be determined from 
the changes in fc£>[24[. Finally, for /i = (or equiva- 
lently kr> = it /(2a)), the characteristic oscillations found 
in this work can be proven rigorously without appealing 
to the perturbation theory [17] ■ In this case, when ig- 
noring the Zeeman coupling, the qp spectrum is gapless 
when kr}£/(2,Tr) is a half-integer and gapped when it is 
an integer. 

Since the effect described here is due to the interference 
effects between the nodal d— wave qps and the ordered 
vortex positions, it should be observable |26j as long as 
the thermal length Lt = hvp/iksT) exceeds I. We ex- 



pect these oscillations to persist in a vortex liquid state 
with strong, albeit only short range, positional order, 
provided that the vortex positional correlation length 
£ p >> Lt » L The true test of these predictions, 
however, would be an observation of the high field quan- 
tum oscillations in the low temperature specific heat in 
the superconducting state with nearly perfect vortex lat- 
tice. The two parameters which control the pattern of 
the oscillations are ku and ao , both of which depend on 
doping and can be determined independently. 
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